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Abstract. A higher order difference equation may be generahy defined in an arbitrary nonempty 
set S as: 

fn(,Xfij Xj-i—lj • • • ) Xn—k) — 9n (^th ^^n— 1 ) • • • ) A;) 

where /„, g^: S^^^^ S are given functions for n = 1, 2, . . . and A; is a positive integer. We 
present conditions that imply the above equation can be factored into an equivalent pair of lower 
order difference equations using possible form symmetries (order- reducing changes of variables). 
These results extend and generalize semiconjugate factorizations of recursive difference equations 
on groups. We apply some of this theory to obtain new factorization results for the important 
class of quadratic difference equations on algebraic fields: 

k k k 

^ ^ ^ ^ 0'i,j,nXn—iXn—j ~\~ ^ ^ bj^nXn—j ~l~ — 0. 
1=0 j=i j=0 

We also discuss the nontrivial issue of the existence of solutions for quadratic equations. 

1 Introduction 

Let 5 be a nonempty set and {fn}'^=i, {gn}'^=i be sequences of functions fn,9n '■ 5*^^^ S where 
k is a positive integer. We call the equation 

fn{Xm ^n— 1 ! ■ ■ ■ i ^n—k^ — 9n{,Xni 3^n— 1 1 • • • j Xn_k) , TL — 1,2,3,... (1) 

a non-recursive difference equation of order k in the set S if fniuo, ■ ■ ■ ,Uk) is not constant in uq 
and gn{uo7 ■ ■ ■ lUk) is not constant in Uk for all n > 1. A (forward) solution of ([T]) is a sequence 
{xn}^=i in S that makes ^ true, given a set of initial values x-j, j = 0, 1, . . . , A; — 1. The recursive, 
non-autonomous equation of order k, i.e., 

Xn — '4'n{,Xji—li • • • ) Xn—k) (2) 
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is a special case of ([T|) with 

fn{uO,Ul, ■ ■ ■ ,Uk) = Uo 

gn{uo,Ui, ... ,Uk) = Ipniui, . . . ,Uk) 

for all n and all uq, ui, . . . , G 5. 

An example of non-recursive equations in familiar settings is the following on the set M of all 
real numbers: 

= a|j;„_i - x„_2|, 0<a<l (3) 

This states that the magnitude of a quantity x„ at time n is a fraction of the difference between 
its values in the two immediately preceding times; however, we cannot determine the sign of Xn 
from (l3|). As a possible physical interpretation of ([3]) imagine a node in a circuit that in every 
second n fires a pulse x„ that may go either to the right (if a;„ > 0) or to the left (if Xn < 0) but the 
amplitude |x„| of the pulse obeys Eq.(l3|). With regard to the variety of solutions for ^ we note 
that the direction of each pulse is entirely unpredictable, regardless of the directions of previous 
pulses emitted by the node; hence a large number of solutions are possible for ([3]). We discuss this 
equation in greater detail in the next section. 

For recursive equations such as ^ on groups, recent studies such as [7], [11], [12], [13], [Tl] , 
|15j . |17j . |18) . show that possible form symmetries (i.e., order-reducing coordinate transformations 
or changes of variables) and the associated semiconjugate relations may be used to break down 
the equation into a triangular system [20] of lower order difference equations whose orders add up 
to the order of ([2]). But in general it is not possible to write ([1]) in the recursive form ([2]) so the 
question arises as to whether the notions of form symmetry and reduction of order can be extended 
to the more general non-recursive context. 

In this paper, we show that for Eq.([TI) basic concepts such as form symmetry and factorization 
into factor and cofactor pairs of equations can still be defined as before, even without a semiconju- 
gate relation. A concept that is similar to semiconjugacy but which does not require the unfolding 
map is sufficient for defining form symmetries and deriving the lower order factor and cofactor 
equations. Using this idea we extend previously established theory of factorization and reduction 
of order to a much larger class of difference equations then previously studied. In particular, we 
apply this extended theory to the important class of quadratic difference equations. 

2 Factorizations of non-recursive equations 

Equation ([1]) generalizes the recursive equation (|2|) in a different direction than the customary 
one, namely, through unfolding the recursive equation of order A; to a special vector map of 
the A;-dimensional state-space S^. Nevertheless, it is convenient to define as states the points 
{xmXn-i, . . . ,Xn-k+i) ^ 5''^ or some invariant subset of it that contains the orbit 
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of every solution {xn}^=i of ([T]). We may designate the point (xq, . . . ,x^k+i) corresponding to 
n = on each orbit as the "initial point" of that orbit; in contrast to recursive equations however, 
there may be many distinct orbits having the same initial point. 

Analyzing the solutions of a non-recursive difference equation such as ([1]) is generally more 
difficult than analyzing the solutions of recursive equations. Unlike the recursive case, even the 
existence of solutions for ([TJ in a particular set is not guaranteed. But studying the form symmetries 
and reduction of order in non-recursive equations is worth the effort. The greater generality of 
these equations not only leads to the resolution of a wider class of problems, but it also provides 
for increased flexibility in handling recursive equations. 

Before beginning the formal study of factorizations of non-recursive equations, let us consider 
an illustrative example that highlights several issues pertaining to such equations. 

Example 1 Let {an}'^=i be any sequence of non-negative real numbers and consider the following 
third- order difference equation on M: 

\Xn ~l~ ^Jfi—il — Oj-ilXji—l Xn— sl' (4) 

By adding and subtracting Xn-2 inside the absolute value on the right hand side of ^ we find that 

\Xn + Xn-l\ = a„|x„_i + Xn-2 " {Xn-2 + a^n-s)! (5) 

The substitution 

tn — Xji -\- Xfi—l (6) 

in (0) results in the second-order difference equation 

\'tn\ = 0.n\tn-l — tn-2\ (7) 

that is related to ^ via ^B^. Eq.^ is analogous to a factor equation for ^ while written as 

Xn — in Xfi—l (8) 

is analogous to a (recursive) cof actor equation. Also the substitution ^ is analogous to an order- 
reducing form symmetry; see, e.g., fWj . fl^ . J7^ or f7g|/ . 

The factor equation ^ is of course not recursive; it is a generalization of in the introduction 
above. If {(3n}'^=i is any fixed but arbitrary binary sequence taking values in {—1,1} then every 
real solution {sn}'^=i of the recursive equation 

Sn = (3nan\Sn-l - Sn-2\ (9) 

is also a solution of This follows upon taking the absolute value to see that {sn}'^=i satisfies 
Eq.^. The single non-recursive equation ^ has as many solutions as can be generated by each 
member of the uncountably infinite class of equations (0j put together. Numerical simulations and 
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other calculations indicate a wide variety of different solutions for (0j with different choices of 
{/3n}^=i and clearly no less is true about Q). 

These facts remain true in the special case mentioned in the Introduction, i.e., equation (3^) in 
which the sequence an = a is constant. By way of comparison first consider the case where = 1 
is constant for all n. Then each solution of the recursive difference equation 

Sn = a\Sn-l - Sn-2|, < a < 1 

is uniquely defined by an initial point (so,s~i) € and 

lim Sn = 0. (10) 

n— ^-oo 

This claim is proved as follows. Without loss of generality assume that sq >0 and define 

fj, = maxjso, si} > 0. 

Then 

S2 = a\si — Sol < amaxjso, si} < a^, 
•S3 = a\s2 — si| < amax{s2, si} < afj,. 

Continuing, 

54 = a\s3 — S2I < amax{s3, 82} < a^/u, 

55 = a|s4 — S3I < amax{s4, S3} < a^fi. 

This reasoning by induction yields 

S2n, S2n+1 < O^fJ; for all U 

which proves [W\) . 

By contrast, if {I3n}'^=i is not a constant sequence then complicated solutions may occur for 
The computer generated diagram in Figure{l\ shows a part of the solution of the difference equation 



with parameter values: 



a/3n|Sra-l — Sn-2\ 



a = 0.8, s__i = So = 1 

-1 ifrn< 0.45 
1 ifrn> 0.45 

3.75r„_i(l - r^-i), ro = 0.4. 
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Figure 1 : A complex solution of equation ([3]) highlighting the unpredictability of sign changes 



The complex behavior seen in Figure{l\is indicative of the stochastic nature of the pulse direction 
in the physical interpretation in the Introduction. Although does not appear explicitly in 

^ the unpredictability of the sign of x„ is a way of interpreting the arbitrary nature of the binary 
sequence. 

Finally, we observe that equation ^ which has order two, also admits a reduction of order as 
follows. If {tn} is a solution of that is never zero for all n then we may divide both sides of ^ 
by \tn-i\ to get 



where the substitution 



1 



tn-2 



tn-1 



t 



n-l 



(analogous to the inversion form symmetry; see JJE/) yields the first- order difference equation 

1 



1 



rn-l 



(11) 



Eq. [Tl\) is related to ^ via the (recursive) equation 



(12) 



The factorizations and corresponding reductions in order given by equations ([7]), ([8]), ([TT]) and 
(I12p are among the types discussed below along with a variety of other possibilities. For additional 
examples and further details see [lOj. 
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2.1 Form symmetries, factors and cofactors 

In this section we define the concepts of order-reducing form symmetry and the associated fac- 
torization for Eq.([T|) on an arbitrary set S. In analogy to semiconjugate factorizations, we seek a 
decomposition of Eq.([T]) into a pair of difference equations of lower orders. A factor equation of 
type 

(t'nitn-i tn—li • • • ■> tn—m) — ''Pniini tn—li • • • i tn—m): 1 ^ TTl ^ k 1 (13) 

may be derived from ([1]) where 4'n,tpn '■ S""^^ — ?■ S for all n if there is a sequence of mappings 
Hn : S'^+i ^ 5"^+! such that 

fn = 4>n° Hn and gn = TpnO Hn (14) 

for all n > 1 . If we denote 

Hn{uo, ... ,Uk) = [/iO,n(uO, • • • , Uk) , hl,n{uo, ■ ■ ■ ,Uk), ... , hm,n{uo, • • • , Ufc)] 

then for each solution of Eq.([T]) 

0n {ho,n (^^n; • • • j Xn—k)i ^l,n (^^n i ■ ■ ■ i ^n—k)i ■ ■ ■ i hm,n (^^nj • • • 5 ^n— fc)) — 
V'n (^0,n (^^nj • • • ; ^n—k^j ^l,ra(^n) • • • ) Xn—k\ • • • ; hm^ni^m • • • ; -^n— fc)) 

In order for a sequence in S defined by the substitution 

tn — h^ nip'^m ■ ■ ■ 1 Xn—k) 

to be a solution of (|13p . the functions must have a special form that is defined next. 

Definition 2 A sequence of functions {Hn} is an order-reducing form symmetry of Eq.mi on 

a nonempty set S if there is an integer m, 1 < m < k, and sequences of functions (j)n,'4'n ■ S"^^^ S 
and hn : 5*^-™+! S such that 

Hn{uo, ... ,Uk) = [hn{uo, ■ ■ ■ ) '^k—m 

), hn-l{ui, . . . 

1 Uk—m+1 Uk)] (15) 

and the sequences {(pn}, {V'n}? {fn}, {dn} o-n-d {Hn} satisfy the relations {l-j^ for all n > 1. 

If S* = (G, *) is a group then Definition [2] generalizes the notion of recursive form symmetry in 
|12j where the components of Hn are defined as 

hn^jiUj, . . .,Uk-m+j) = Uj * hn^j{Uj+i, . . . ,Uk-m+j), j = 0,1, . . . ,m. 

We have the following basic factorization theorem for non-recursive difference equations. 
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Theorem 3 Assume that Eq.^^ has an order-reducing form symmetry {Hn} defined by il5\} on a 
nonempty set S. Then the difference equation ^^ has a factorization into an equivalent system of 
factor and cofactor equations 

4'n{tm • • • 5 tn—m) — V'n(^rt) • • • ) ^n— m) (-'-^) 
hn{Xii^ • • • 1 Xn—k+m) — (1'^) 

whose orders m and k — m respectively, add up to the order of (CP- 

Proof. To show the equivalence, we show that for each solution of ([T]) there is a solution 

{(*ni2/n)} of the system of equations PT|) such that y„ = Xn for all n > 1 and conversely, for 
each solution {{tn,yn)} of the system of equations (fT6|) and pT|) the sequence {yn} is a solution of 
©• 

First assume that {xn} is a solution of Eq.([T|) through a given initial point (xq, X-i, . . . , x_fc+i) G 
(^fc+i^ Define the sequence {tn} in S" as in (117p for n > — m + 1 so that by (I14p 



</'n(^n) • • • ) tn—m) — 4'n{f^n{^n} • • • 5 . . . , hn—m{p^n—mi ■ ■ ■ 1 ^n—k)) 

— (pn^Hni^Xn^ • • • ; ^n— A;)) 

— fni^m • • • , ^^n— fc) 

— 9n{Xm • • • 1 Xn—k) 

— '4'n(,Hn{Xm • • • , S^n— fc)) 

— "^nitni ■ ■ ■ 1 tn—m)- 

It follows that is a solution of (|16p . Further, if y„ = x„ for n > —k + m then by the 
definition of t„, {yn} is a solution of (fT7|l . 

Conversely, let be a solution of the factor-cofactor system (fT6]l . (fT7|l with initial values 

ioi • • • > i-m+ii y-m, • • • y-fc+i ^ C. 

We note that yo, . . . , y^m+i satisfy the equations 

hjiyj, . . .,yj-k+m) = tj, j = 0,-1,. . . , -m + 1. 

Now for n > 1, (I14p implies 

/n(yn, • • • ,yn-k) = 4>n{Hn{yn, ■ ■ ■ .Vn-k)) 

— 't'nihniyni ■ ■ ■ 1 Un-k+m)-! ■ ■ ■ 1 hn—miyn—mt • • • 1 Un-k)) 

— 't'n (tji , • • • , in— m) 

— V'n(in; • • • ) tji—m) 

— ''Pn{Hn{Xni • • • ; ^n-k)) 

— dni^ni • • • 1 X-a-k)- 
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Therefore, {un} is a solution of ([T]). ■ 

The concept of order reduction types for non-recursive difference equations can be defined sim- 
ilarly to recursive equations as in prior studies and is not repeated here. 

2.2 Semi-invertible map criterion 

A group structure is necessary for obtaining certain results such as an extension of the useful 
invertible map criterion in [12] and [15] to non-recursive equations. In this section we assume that 
S = {G, *) is a non-trivial group with the goal of obtaining an extension of the invertible map 
criterion to the non-recursive equation (jlSp . For a discussion of some issues pertaining to difference 
and differential equations on algebraic rings see [1]. 

Denoting the identity of G by l, the difference equation ([T|) can be written in the equivalent 
form 

En{Xn, Xn^i, X^-k) = h > 2, n = 1, 2, 3, . . . (18) 

where En = fn * [9n]^^ (the brackets indicate group inversion). A type-(/c — 1, 1) reduction of 
Eq. ()18p is characterized by the following factorization 

4>n{tn, ■ ■ ■ ,tn-k+l) = 1^ (19) 
hn{Xn,Xn-l) =tn (20) 

in which the cofactor equation has order one. This system occurs if the following is a form symmetry 
of (dHD: 

Hniuo, ...,Uk) = [hn{uo,Ui),hn-l{ui,U2), ■ ■ ■ , /l„_fc+l (ufc-1 , ""A:)] • (21) 

A specific example of the above factorization is the system of equations ([7]) and ([8]) in Example 
[1] which factor Eq.(j3|). More examples are discussed below. 

The following type of coordinate function hn is of particular interest in this section. 

Definition 4 A coordinate function h : G'^ ^ G on a non-trivial group G is separable if 

h{u, v) = /i(n) * 9{v) 

for given self-maps fi,9 of G into itself. A separable h is right semi-invertible if 6 is a bijection 
and left semi-invertible if ^ is a bijection. If both 6 and fi are bisections then h is semi- 
invertible. A form symmetry {Hn} is (right, left) semi-invertible if the coordinate function hn is 
(right, left) semi-invertible for every n. 

Note that a semi-invertible h is not a bijection in general; for instance, consider h{u, v) = u — v 
where G is the group of all real numbers under addition. 

Clearly, functions of type u*h{v) where /i is a bijection are semi-invertible functions. Therefore, 
semi-invertible functions generalize the types of maps discussed previously in [12], [13], [E] and 
|18| . The next theorem shows that the invertible map criterion discussed in these references can be 
extended to all right semi-invertible form symmetries. 



8 



Theorem 5 (Semi-invertible map criterion) Assume that hn{u,v) = /in('u) * On{v) is a sequence 
of right semi-invertible functions with bijections On of a group G. For arbitrary uq,vi, ■ ■ ■ ,Vk G G 
define Co,n = uq and for j = 1, . . . , k 

Cj,ni^O,Vl, ■■■,Vj) = 6'^_^j+i([/in-i+l(Ci-l,n(^iO,^l,- • * Vj) (22) 

with the usual distinction observed between map inversion and group inversion. Then Eg. U8\) has 
the form symmetry i21\) and the associated factorization into equations il9\) and l\20\) if and only 
if the following quantity 

Eniuo,Cl,ni'>^0,Vi), . . . , Ck,ni'^0, Vl, ■ ■ • (23) 

is independent of uq for all n. In this case, the factor functions (f>n 

are given by 

0„(V1, ...,Vk)= Eniuo,Cl,n{'^0,Vl), ■ ■ ■ , Ck,ni'^0, Vl, ■ ■ -^Vk))- (24) 

Proof. Assume that the quantity in ()23p is independent of uq for every n. Then the functions (pn 
in (p^ are weU defined and if Hn is given by (j2T]) and Wj+i = hn^j{uj, Uj+i) for j = 0, . . . , A: — 1 in 
(I22D then 

4>n{Hn{uO, ■ ■ ■ ,Uk)) = 4'n{hn{uo,Ui),hn-l{ui,U2), ■ ■ ■ , hn-k+l{uk-l,Uk)) 
= En{uo, Cl,n{uO, hn{uo, Ui)), ... , 

C,k,n{uQ,hn{uQ,Ui),. . . ,hn-k+l{Uk-l,Uk)). 

Now, observe that 

Ci,n{'^o,Kiuo,ui)) = en^ilnniuo)]''^ * /i„(no) * 6'„(ni)) = Ui. 
By way of induction, assume that for j < k 

Cj,n{uO,hn{uO,Ul), . . . ,hn-j+l{Uj-i,Uj)) = Uj (25) 

and note that 

0+l,n(wO,- • . ,hn-j{Uj,Uj+i)) = e~^-{[^ln-j{Cj,n{uo, ■ ■ . , hn-j+l{Uj-i, Uj)))]^'^ 

* lln-j{Uj) * 9n-j{uj+l)) 
= K-ji\l^n-j{Uj)]'^ * fJ-n-jiuj) * On-jiUj+i)) 
= Uj+i. 

It fohows that (p5]) is true for ah j = 0, 1, . . . , and thus 

(t>n{Hn{uo, . . . , Uk)) = En{uo, ■ ■ ■ ,Uk) 
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i.e., {Hn} as defined by by (j2ip is a form symmetry of Eq. (jl8p and therefore, Theorem [3] implies 
the existence of the associated factorization into equations ([T9|) and (pO]) . 

Conversely, suppose that {Hn} as defined by by (j21|) is a form symmetry of Eq. (jl8p . Then there 
are functions (pn such that for all uq,vi, . . . ,Vk £ G 

En{uo, Cl,n, • • • , Ck,n) = 4>n{Hn{uo., Cl,n, • • • , Cfc.n)) 

= (pnihniuo, Cl,n), ■ ■ ■ , ^n-fc+l (Cfc-l.n, Ck,n)) 

where Cj,n = CjA'^o, vi, . . . , vj) for j = 1, . . . , k. Since 

hn-j+l{Cj~l,n, Cj,n) = Mn-j+l (Cj-l.n) * ^n-j+l(Ci,n) 

= Mn-i+l(Cj-l,n) * 6l„_j + i(6l~^_^i([/i„_j+i(Cj-l,„,)]""^ * Vj)) 
= fJ-n-j+l{Cj-l,n) * [^J'n-j+l{Cj-l,n)]~'^ * Vj 



it follows that £'„(uo, Ci.n; • • • ) Cfc.n) is independent of uq for all n, as stated. ■ 

Special choices of /x„ and give analogs of the identity, inversion and linear form symmetry 
that are discussed in [18]. The next example illustrates both Theorem [5] and the significant fact 
that a recursive difference equation may have non-recursive form symmetries. 

Example 6 Consider the recursive difference equation 



aXn_i + bXn-l + CXn-2 + d (26) 

where a, b,c,d(z R. To find potential form symmetries of this equation, first we note that every real 
solution of i26\) is a (non-negative) solution of the following quadratic difference equation 

x\ - O-X^-i - bXn-l - CXn-2 - d = 0. (27) 



Based on the existing terms in (21), we explore the existence of a right semi-invertible form sym- 
metry of type 

h{u,v) = + av. (28) 

We emphasize that (28\) is not a recursive form symmetry of the type discussed in (TBI or fTSj. Here 

9^^{v) = v/a so 

Ci(«o,^i) = -(-^^o + ^i), (29) 
a 

C2(no,ui,f2) = --"2 t(-Wo + ui)^. (30) 
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By Theorem 128\) is a form symmetry for i26\) if and only if the following quantity is independent 
ofuQ 

Eiuo, (1X2) =4- aCf - bCi - 0(2 - d. (31) 

Using i29\) and \30\) in i31\) and setting the coefficients of all terms containing uq equal to zero 
gives the following two distinct conditions on parameters 

b a c 

l + - = 0, - + — =0. 

From these we obtain 

a = —b, c = aa = —ah. (32) 

If b,c ^ then conditions [32\) imply the existence of a form symmetry — bv for [2^) with a 
corresponding factorization: 

tn - atn-1 = d, (33) 

xl - bXn-1 = tn- (34) 

The positive square root of x„ in the cof actor equation |3^[ ) can be used to obtain a factorization of 
the recursive equation i26\) as 

tn = atn-i + d, 

Xn = \JbXn-\ + in, ^0 = 3:0- bX-X- 

Note that this factorization of l^26\) as a system of recursive equations is derived indirectly via 
the non-recursive equation |g7[ ) and its factorization 133\) and l^34\ ) rather than directly through a 
semiconjugate relation. 

The existence and asymptotic behaviors of real solutions discussed in the preceding example 
are not as easily inferred from a direct investigation of (j26p . For a discussion of solutions of (j26p 
using the above factorization see [lOj. 

If {Hn} is a semi-invertible (right and left) form symmetry of Eq. ljlSp then the following result 
states that the cofactor equation (|20|) can be expressed in recursive form. 

Corollary 7 Assume that the functions hn in Theorem\^are semi-invertible so that both fin CLnd 
9n are bijections. Then Eq. has the following factorization 

(pnitn, ■ ■ ■ , tn-k+l) = 
Xn = fJ-n^itn * [dniXn~l)]~^) 

in which the cofactor equation is recursive. 
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3 Quadratic difference equations 

In [12j it is shown that a non-homogeneous and non-autonomous hnear difference equation 

Xn + ai^nXn-1 H h ak,nXn-k = K, ak,n 7^ for all n > (35) 

has the linear form symmetry and with the corresponding SC factorization over a non-trivial alge- 
braic field T if the associated Riccati equation 

Oin = aO,n H \ \ h 



a„_l an-\an-2 C^n-l ■ ■ ■ Otn-k 



of order k — 1 has a solution {a„} in It can be checked that this condition is equivalent to the 
existence of a solution of the homogeneous part of ([35]) 



Xn + ai^nXn-1 H h ttk^nXn-k = (36) 

that is never zero. For if {y™} is a nonzero solution of ([36]) then the ratio sequence {yn+i/Un} is a 
solution of the Riccati equation above. These facts lead to a complete analysis of the factorization 
of (|35]) over algebraic fields; see [10] for additional details. Riccati difference equations have been 
studied in [9] (order one) and [3] (order two). 

A natural generalization of the linear equation (|35p is the quadratic difference equation over a 
field ^ 

k k k 
^ ^ ^ ^ 0'i,j,nXn—iXn—j ~l~ ^ ^ bj^nXn—j ~l" — (37) 

that is defined by the general quadratic expression 

k k k 

CI'i,j,nUiUj + y ^ bj^nUj + Cn- 
1=0 j=i j=0 

Linear equations are obviously special cases of quadratic ones where aij^n = for all i,j,n. 
Further, equations of type (j37p also include the familiar rational recursive equations of type 

Ek ^^k ^^k 1 

_ 1=1 l^j=i 0'i,j,nXn-iXn-j — Z^j=l Oj,nXn-j — Cn 

2^j=l 0'(},j,nXn—j 1 OQ^n 

as special cases where 

oo,o,n = for all n. (39) 
Special cases of (I38p over the field M of real numbers include the Ladas rational difference 



equations 

^ + J2j=l ^j^n-j 
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as well as familiar quadratic polynomial equations such as the logistic equation 



the logistic equation with a delay, e.g. 



Xn = Xn-l{a - bXn-2 




and the Henon difference equation 



a + bx, 




n— 1- 



These rational and polynomial equations have been studied extensively; see, e.g., [2], [3], [5], 



We note that all solutions of ([38j) are also solutions of the quadratic ([37|) when condition ([39]) 
is satisfied. The extensive and still far-from-complete work on rational equations of type (j38p is a 
clear indication that unlike the linear case, the existence of a factorization for Eq. ()37p is not assured 
and in general, finding any factorization into lower order equations is a challenging problem. 

3.1 Existence of solutions 

When (j39p holds solutions of (j38p may be recursively generated. Difficulties arise only when the 
denominator becomes zero at some iteration; such singularities often arise from small sets in the 
state-space and are discussed in the literature on rational recursive equations. However, if 



then the problem of the existence of solutions for (j37p is entirely different in nature. If (|40p holds 
then the issue is not division by zero but the existence of square roots in the field T. We examine 
this issue for the familiar field M of real numbers. The next example illustrates a key idea. 

Example 8 Let a, 6, c be real numbers such that 



m, m, m, m- 



oo,o,n 7^ for all n 



(40) 



a ^ and b, c 



> 0. 



(41) 



and consider the quadratic difference equation 



(42) 



We solve for the term x„ by completing the squares: 



z I z 



2 I A 2 , 
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Now we take the square root which introduces a binary sequence {j3n}'^=i with j3n S {—1, 1} chosen 
arbitrarily for every n: 



Xn = -Xn-l+ Pn\l^xl_^+bxl_^ + C (43) 

Under conditions (^Tp, for each fixed sequence {j3n}'^=i every solution of the recursive equation (^^p 
with real initial values is real because the quantity under the square root is always non-negative. 
Furthermore, since 



^Xn-l 



it follows that for each n, 



Xn> if I3n = 1, 
x„ < if I3n = -1. 

This sign- switching implies that a significant variety of oscillating solutions are possible for Eq. 
under conditions (fTip - Indeed, since f3n is chosen arbitrarily, for every sequence of positive integers 

{mi, 7712, ms, • • ■} 

there is a solution of ^^'^^ starts with positive values of Xn for mi terms by setting /3„ = 1 for 
1 < 71 < mi. Then x„ < for the next m2 terms with /3„ = — 1 for n in the range 

mi + 1 < 77 < 7771 + "72 

a77d so on with the sign of x^ switching according to the sequence {mn}'^=i. 
Figure [2] illustrates the last part of the above example. 

The method of completing the square discussed in Example [8] can be applied to every quadratic 
difference equation with real coefficients. This useful feature enables the calculation of solutions of 
such quadratic equations through iteration, a feature that is not shared by non-recursive difference 
equations in general. The next result sets the stage by providing an essential ingredient for the 
existence theorem. 

Lemma 9 /77 the quadratic difference equation 

^ ^ ^ ^ OiiJ^nXn—iXn—j ~l~ ^ ^ bj,nXn—j ~l~ — (44) 
i=0 j=i j=0 
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Figure 2: Oscillations in each solution x„ exhibit the same pattern as those in the binary sequence 
generating it 



assume that all the coefficients aij^n,bj^n,Cn are real numbers with ao,o,n 7^ for all n. Then 

is a real solution of ( (^^[ j if and only if is a real solution of some member of the family of 

recursive equations, namely, the recursive class of 

) + Qni ) (45) 

where {/3„} is a fixed hut arbitrarily chosen binary sequence with values in the set { — 1,1} and for 
each n, 

L„(ui, ...,Uk) = ^ ( 6o,n + (46) 



-I ( ' \ 

),0,n \^ ^.^^ J 

''0'" \i=l j=i j=l 



Qn{ui,..., Uk) = 17,7, aij^nUiUj + } ^ bj-nUj + Cn \ ■ (47) 

«0,0,r 



Proof. If 



k k k 
En{uo,Ui,. • • , ""fc) = ^ ^ ai,j,nUiUj + ^ bj^nUj + Cn 
1=0 j=i j=l 

then using definitions (j46p and (j47p we may write 

En{uo,Ui, ...,Uk)= ao,o,nkn " '^Ln{ui, . . .,Uk)Xn - Qn{ui, ■ ■ • , Wfc)] (48) 

Since ao,o,n 7^ for all n, the solution set of (|44p is identical to the solution set of 

- 2Ln{ui, . . .,Uk)Xn - Qn{ui, . . . ,Ufc) = 0. (49) 
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Completing the square in (jl9]) gives 

[Xn - Lniui, . . .,Uk)f = Ll{ui, . . . , Ufc) + Q„(ui, . . . 

which is equivalent to (j45p . ■ 

Example [8] provides a quick illustration of the preceding result with 

L{ui,U2) = ^ui, Q{ui,U2) = bul + c 

where Ln = L and Qn = Q are independent of n. 

We are now ready to present the existence theorem for real solutions of (l44|l . Let {Z?™} be a 
fixed but arbitrarily chosen binary sequence in {—1,1} and define the functions L„ and Qn as in 
Lemma [9l Further, denote the functions on the right hand side of (j45p by fn, i.e., 

fn{ui, ...,Uk) = Ln{ui, . . . ,Uk) + (in\/ Ln{ui, . . . , Ufc) + Q„(til, . . . ,Uk). (50) 

These functions on M'^ unfold to the self-maps 

F„(ui, ... ,Uk) = {fn{ui, . . . ,Uk),Ul, . . . ,Uk~l) ■ 

We emphasize that each function sequence {Fn} is determined by a given or fixed binary se- 
quence {/3n} as well as the function sequences and {Qn} that are given by (j44|l . Clearly the 
functions fn are real- valued at a point (ui, . . . , Uk) G if and only if 

L^(ui, . . . , Ufc) + Qn(^^i, . . . , -Ufc) > 0. (51) 

Theorem 10 Assume that the following set is nonemtpy: 

oo 

5= Pi {{ui,...,Uk) eR'' : Ll{ui,...,Uk) + Qn{uu...,Uk)>oY (52) 

n=0 

(a) The quadratic difference equation ftas a reaZ solution {xn}'^=_k+i '^'^^ '^'^^2/ ^/ poini 
-Pq = (s^O) 2:^-1) •••) 2;_^.+i) is in S and there is a binary sequence {/3n} in {~^A} such that the 
forward orbit of Pq under the associated maps {Fn} is contained in S; i.e., 

{F„oF„_io...oFi(Po)}^=i C5. 

(b) If the maps {Fn} have a nonempty invariant set M d S for all n, i.e., 

Fn{M) CM CS 
then the quadratic difference equation (44\ ) has real solutions. 
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Proof, (a) By Lemma [9l is a real solution of (j44p if and only if there is a binary 

sequence {/3n} in {—1,1} such that {xn}^^_i^_^_i is a real solution of the recursive equation ()15]) . 
Using the notation in (j5U|) . Eq. (j45|) can be written as 

Sra = fn{Sn—l: • • • ? ^n—k)- (^3) 

Now the forward orbit of the solution of ([53|) is the sequence 

O = {F„ o F„_i o . . . o Fi(Po)}~ 1 = {(X„-1, . . . , Xn-k)}n=l 

in R*^ that starts from Pq E 5. It is clear from the definition of T that each x„ is real if and 
only if O C <S. This observation completes the proof of (a). 

(b) Let Po e M. Then {F„ o Fn-i o • • • o Fi(Po)}~=i C M C 5 so by (a) 

{{Xn-li • • • ) 2;n-fc)}J^Ll C S. 

It follows that each Xn is real and thus, of ([53]) . An application of Lemma [9] now 

completes the proof. ■ 

If S has any invariant subset M = M({/3„}) (relative to some binary sequence {/3n} in {— 1) 1} 
or equivalently, to some map sequence {Fn}) then the union of all such invariant sets in S is again 
invariant relative to all relevant binary sequences {/?„} (or map sequences {Fn}). Invariant sets 
may exist (i.e., they are nonempty) for some binary sequences {/3n} in 1} and not others. 
However, the union of all invariant sets, 

M = \^{M{{/3n}) ■ {/3n} is a binary sequence in {—1, 1}} 

is the largest or maximal invariant set in S and as such, A4 is unique. In particular, if S is invariant 
relative to some binary sequence then A4 = S. 

We refer to A4 as the state-space of real solutions of (j44p . The existence of a (nonempty) A4 
may signal the occurrence of a variety of solutions for ()44p . In Example El where 5 = R'^ is trivially 
invariant (so that A4 =B?) we observed the occurrence of a wide variety of oscillatory behaviors. 
Generally, when 5 = R'^ every solution of (I44p with its initial point Pq € R'^ is a real solution. The 
next result presents sufficient conditions that imply 5 = R'^. For cases where 5 is a proper subset 
of R'= see 

Corollary 11 The state-space of real solutions of |^^[ ) is R*^ if the following conditions hold for 
all n: 

Aj^n > 0, Bij,^ = 0, (54) 

k fj2 

^ -j^ < bi n - 4ao,o,nCn (55) 
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where for k > 2, j = 1, . . . ,k, and n > 

-Bi,j,n — O,0,i,nXl0,j,n 2fflo, 0, 71*^2, jr,n) ^ ^ jj 
Cj^n = aOj,nbo,n — 2ao,0,nfej,n- 

Proof. By straightforward calculation the inequality L'^ + Qn > 0, i.e., ()5ip is seen to be equivalent 
to 

k k—1 k k 

+ ^ X] X] Bi,j,nUiUj + 2 ^ C'j,nlij + &0,n > 4ao,0,nCn (56) 

i=i «=i i=«+i j=i 

for all (ui, . . . , Uk) G with the coefficients Aj^n, Bij^n and Cj^n as defined in the statement of the 
corollary. By conditions the double summation term in drops out and we may complete 
the squares in the remaining terms to obtain the inequality 



k y Q \2 (7? 



By ()54p the left hand side of the above inequality is non-negative while its right hand side is 
non-positive so (j56p holds under conditions (j54p and (j55p . The proof is completed by applying 
Theorem [lOl ■ 

The next result, which generalizes Example [8l is obtained by an immediate application of the 
above corollary to the non-homogeneous quadratic equation of order two with constant coefficients. 

Corollary 12 The quadratic difference equation 

xl + ao,lX„Xn-l + ao^2XnXn-2 + ai.l^^n-l + ai,2Xn-lXn-2 + 
+a2,2xl_2 + boXn + biXn-1 + b2Xn-2 + C„ = 

has as a state-space of real solutions if the following conditions are satisfied: 

^ "-0,1 ^0,2 1 J 

ai,i < ^p, 02,2 < — , ai,2 = 2"o,iao,2 and 



^ bl {ao,ibo - 26i)2 (ao,2?>o - 262)^ . 
Cn< -V - -rr^ r - -77-^ : r for all n. 

" - 4 4(a2^ -4ai,i) A{al^-4a2,2) 



3.2 Factorization of quadratic equations 

In the remainder of this paper we investigate conditions for the possible existence of a special semi- 
invertible form symmetry for quadratic difference equations in which both fin and On are linear 
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maps on T for all n; hence, this is called a linear form symmetry. This is a first step in a broader 
study of the latter type of equation on algebraic fields; see [10] for further discussion. 

To simplify the discussion without losing sight of essential ideas, we limit attention to the case 
k = 2, i.e., the second-order equation 

En{Xn,Xn-l,Xn_2) = (57) 

on a field J- where 

En{uo,Ui,U2) = ao,o,n^to + ao,l,nUoUi + ao^2,nUoU2+ (58) 
ai,0,nUl + ai,2,nUlU2 + a2,0,nU2 + 
bo,nUO + bl,nUl + &2,n^i2 + C„. 

In this expression, to assure that uq (corresponding to the Xn term) does not drop out, we may 
assume that for each n, 

if ao,o,n = ao,i,™ = ^o,n = or ai,2,n = a2,o,n = b2,n = 0, then ao,2,n / 0. 
If the quadratic expression (j58p has a form symmetry {Hn} where 

Hn{uo,Ui,U2) = [hn{uo,Ui),hn^l{ui,U2)] 

then there are sequences {(pn}, {hn} of functions (/)„, hn : T"^ ^ T such that for all (uq, ni, 1^2) G T^, 

(t>n{hn{uo,Ui),hn-l{ui,U2)) = £'„ (uo, Ui , U2) . (59) 

The form symmetry {Hn} is semi-invertible if there are bijection ^^(^n '■ T ^ T such that 

hn{u,v) = fln{u) + On{v). 

The next corollary of Theorem [5] presents conditions that imply the existence of a linear form 
symmetry for Eg. ([5^ 



Corollary 13 The quadratic difference equation ( |57[ ) has the linear form symmetry 

hn{u, v) = OnU + V, a„ 7^ for all n (60) 

if and only if a sequence {an} exists in the field T such that all four of the following first-order 
equations are satisfied: 

2 2 2 2 / \ 

aO,0,ra — aCl.raOn + Q-\,\,nOLn + a0,2,n«nara-l " ai,2,n«nC^ra-l + 0'2,2,nOLn(^n-\ =0 (^1) 

aO,l,n - 2a\XnOin - a0,2,nan-l + 2a\;2,nOin(y^n-\ - 2a2^2,nOinC?n-\ = (62) 

ao,2,n - ai,2,nan + 2a2,2,nO:nan-l = (63) 

bo,n - bi^nOtn + &2,nanan-l = 0. (64) 
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In this case, Eg. |57| j has a factorization with a first- order factor equation 

(ai,l,n - ai,2,nan-l + a2,2,na„_l)in + «2,2,nin-l + (ai,2,n " 2a2,2 
+ (&l,n - &2,nan-l)in + &2,nin-l + C„ = 0. 

and a cofactor equation anXn + a^n-i = o^feo 0/ order one. 

Proof. For any nonzero sequence {«„} in T the functions /i„ defined by (j60p are semi-invertible 
witli ^„(u) = a^u and ^nl'i') = "^^ for all n. Note that = On for all n and the group structure is 
the additive group of the field so the quantities Cj> in Theorem [5] take the forms 

Cl,n = Cl,n{'^0,Vi) = -a„Mo + ^'l, 

C2,n = C2,n{uO,Vl,V2) = a„a„-lUo - On-lfl + V2. 

By Theorem [Sj Eq. (j57p has the linear form symmetry if and only if the expression En = 
Eniuoi Ci,n) C2,n) is independent of uq for all n. Now 

En = ao,o,nWo + aQ,i,nUo{vi - anUo) + 

aQ^2,nUo{anan~lUQ - an^lVl + V2) + ai^nivi - anUQf + 
ai,2,n{vi - anUo){anan-lUQ - Un-lVi + V2) + 

'^2,2,71 

bi,n{vi - a„no) + h2,n[anan-lUo - an-\V\ + V2) + Cn- 
Multiplying terms in the above expression gives 

En = aOfl,nul + OQ^i^nUQVi - ao,l,nanMo + «0,2,nanan-lMo " 
a0,2,nan-lU0Vi + aQ^2,nU0V2 + tti^i^nvl - 2ai,i,„a„uo't^i + 
ai,l,na^n'^l + a'l,2,nanan^lU0Vi - ai^2,nan-ivl + ai^2,nVlV2- 
ai,2,nO?n(^n-lul + ai^2,nanan-lUoVi - ai^2,nanUoV2 + 
a2,2,nOlnO?n-iul - 2a2,2,nanan_l ^O^^l + 2o2,2,nanan-l'"0^'2 + 

{an-lVl - V2f + 6o,nMO + &l,nf 1 " ^l,nanMO + 
h2,nanan-lU0 - 62,n(an-l^^l " V2) + Cn- 

Terms containing uq or Mq must sum to zeros. Rearranging terms in the preceding expression 
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gives 



En = (aO,0,n — aO.l.nOn + ao,2,nan«n-l + Ql.l.nOn ~ 0'\,2,nOLn'^n-\^ 

a2,2,na^an_i)'Uo + ('^0,l,n — ao,2,n«n-i " 2ai^i^na„+ 
2ai,2,nanan-i - 2a2,2,nanan-i)^owi + (ao,2,n - ai,2,nan+ 

2a2,2,nanan-l)^i0^^2 + (&0,n " ^x^nOin + &2,nanan-l)'Uo + 
(ai,l,n - ai,2,nan-l + a2,2,nan-l)vj + a2,2,nv\ + (ai,2,n- 
2a2,2,nan-l)'yi'y2 + (^l.n " fe2,nan-l)^^l + b2,n^^2 + Cn- 

Setting the coefficients of variable terms containing uq equal to zeros gives the four first-order 
equations ([6T]) - ([6l|) . The part of En above that does not vanish yields the factor functions 

(l>n{vi,V2) = (ai,l,n - ai,2,nan-l + a2,2,na\_i)vl + a2,2,n^2 + 

(ai,2,n - 2a2,2,nan-l)?^l?J2 + (&l,n " &2,nan-l)'yi + 62,n'V2 + Cn- 

This expression plus the linear cofactor tn = hn{xn, Xn-i) = otnXn + Xn-i give the stated 
factorization. ■ 

An immediate consequence of Corollary 1131 is that every second-order, non- homogeneous linear 
difference equation 

Xn + bl^nXn~l + b2,nXn-2 + C„ = 

has a linear form symmetry and the corresponding factorization into a pair of equations of order 
one (also non- homogeneous, linear) if and only if the first-order difference equation (j64p has a 
solution {an} in J-'\{0}. The existence of a linear form symmetry for non-homogeneous linear 
equations of all orders can be established by a calculation similar to that in the proof of Corollary 
[T3l However, as noted earlier, a complete proof for the general case is already given in jl2] using the 
semiconjugate factorization method which applies to linear equations because they are recursive. 
Therefore, we need not consider the linear case any further here. 

If all coefficients in En{uo,ui,U2) are constants except possibly the free term Cn then a simpler 
version of Corollary 1131 is obtained as follows. 

Corollary 14 The quadratic difference equation with constant coefficients 

aOfixl + ao,lXnXn-l + ao,2XnXn-2 + ai,i3;^_i + ai,2X„_iX„_2 

+02,22:^-2 + boXn + biXn~l + b2Xn-2 + C„ = (65) 

in a non-trivial field T has the linear form symmetry with h{u, v) = au + v if and only if the 
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following polynomials have a common nonzero root a in T : 

00,0 — 00,1" + (oi.i + ao,2)«^ — 01,2"^ + a2,2«^ = 0, (66) 

oo,i - (ao,2 + 2ai^i)a + 2ai,2a^ - 2a2,2«^ = 0, (67) 

ao,2 — ai,2a + 202,20^ = 0, (68) 

60 - hxa + h2C? = 0. (69) 

// such a root a ^ exists then Eq. \63^) has the factorization 

(ai,i - ai,2a + a2,2a^)in + (01,2 - 2a2,2a)intn-i + 02,2*^-1 
+{hi - b2a)tn + h2tn~l + = 0, 

_ _ 1 ,tn 

a a 

Equalities ()66p -()69p often lead to suitable parameter restrictions implying the existence of a 
linear form symmetry for a given difference equation. Here is a sample. 

Example 15 Consider the following quadratic equation (no linear terms) 

lXn~2 = CTn- (70) 



In the absence of linear terms in |7Q[ ) equality i69\) holds trivially; the other three equalities /i66]) - A68]) 
take the following forms 

1-aa + ba^ - ca^ = (71) 
a-ba + 2ca^ = (72) 
6 - ca = 0. (73) 



From (73) it follows that a = b/c. This nonzero value of a must satisfy the other two equations in 
the above system so from \7S^ we obtain 

a + — = 

c 

while (fTiy yields 

ab 

1 = or c = ab. 

c 

Eliminating b and c from the last two equations gives 

b = — a^, c = —a^ and a = —. 

a 
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These calculations indicate that the quadratic equation 

2 2 3 

-\- aXfiXji—X CI XjiXfi—2 Oj Xn— 1X^—2 — (J-ii 

has the linear form symmetry with the corresponding factorization 

a t^ ^ tntn—l — '^ni 
Xn ^Xyi — \ — atyi . 

Corollary (Tl] also implies the non-existence of a linear form symmetry. The next example 
illustrates this fact. 

Example 16 Let a, 6 € C with b ^ and let {c„} be a sequence of complex numbers. Then the 
difference equation 

xi + axl_i + bxl_2 = Cn (74) 

does not have a linear form symmetry because the equality /i6S^) in Corollary\T^does not hold for a 
nonzero complex number a. On the other hand, the substitution yn = x"^ transforms \74\) into the 
non-homogeneous linear equation 

Un + ay-n-l + by n-2 = Cn 

which does have a linear form symmetry in C (by Corollary \14\ with equality ( 1^9)) implying that —a 
is an eigenvalue of the homogeneous part). Substituting for yn in the resulting cofactor equation 
gives a factorization of |7-^[ ). 
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